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Abstract 

We study the influence of a possible H dibaryon condensate on the equation 
of state and the overall properties of neutron stars whose population otherwise 
contains nucleons and hyperons. In particular, we are interested in the question 
of whether neutron stars and their masses can be used to say anything about 
the existence and properties of the H dibaryon. We find that the equation of 
state is softened by the appearance of a dibaryon condensate and can result in 
a mass plateau for neutron stars. If the limiting neutron star mass is about 
that of the Hulse- Taylor pulsar a condensate of H dibaryons of vacuum mass 
~ 2.2 GeV and a moderately attractive potential in the medium could not be 
ruled out. On the other hand, if the medium potential were even moderately 
repulsive, the H, would not likely exist in neutron stars. If neutron stars of 
mass ~ 1.6Mq were known to exist, attractive medium effects for the H could 
be ruled out. 



1 Introduction 



Since Jaffe proposed that there may exist a stable dihyperon (a quark composite with 
baryon number two) Q, an ongoing quest for this particle began [§]. Recent searches 
using kaon beams or heavy ion beams ||, ||, |6| found no candidates or are still 
in progress f7j. There exist some claims for evidence for the H dibaryon produced 
in proton-nucleus || and in heavy-ion collisions ||. Nevertheless, these candidates 
might be misidentified K\ as seen in |J. For a most recent overview on the search 
for the H dibaryon we refer to 



10] 



There are numerous mass estimates for the H dibaryon and they are reviewed in 
IT. | . The existence or nonexistence of the H dibaryon is strongly connected with the 



observation of double A hypernuclei which has been discussed in [E| . Three double A 
hypernuclei have been reported in literature: 



AA 



He 



h AA 



10 Be |4 



and B 



15, 16 



The two A's can decay by strong interactions to the H dibaryon. As this has not been 
seen in the above events, the H must either be heavier than niu > 2m\ + 5aa ~ 2.22 
GeV [117] or the events are misidentified as an H hypernucleus with a shallow attractive 
nuclear potential |lT |. A more stringent condition is the observation of the weak 

2190 
-3.1 



mesonic decay of the double A hypernuclei giving m# > m^ + rrip + m^- + 



MeV [[L8|] where B\ depends on the mass of the decay fragment and is Bp 



1 



MeV for |He. In all deeply bound H dibaryon seems to be ruled out by these 

events. 

If the H dibaryon exists, it will have a certain impact also on the properties 
of dense matter. It is quite established nowadays, that neutron stars have a large 
hyperon fraction in the core and might be described as giant hypernuclei, though 
bound by gravity [fjj. Here again, the presence of hyperons might restrict certain 



properties of the H dibaryon. Recently, studies for neutron stars have been done 
for nuclear matter without hyperons but including H dibaryon condensation |2(J and 
limits have been set for the coupling constants of the H dibaryon |f2T| . 

There might exist heavier partners of the H dibaryon, lumps of strange quark 
matter dubbed strangelets. There are several heavy-ion experiments dedicated to 
search for this novel form of matter |j2"2 , |23| , 23]. In the MIT bag model, strangelets 



with A < 6 are found to be unbound [25|. Nevertheless, light strangelet candidates 



in the range of 6 < A < 40 might be stable against weak hadronic decay |26|, |27J (an 



overview of the properties of strange matter can be found in [28]). The H dibaryon as 



well as these light strangelets can occur in dense matter as a precursor of the phase 
transition to a quark plasma. 

In this paper, we study the influence of H dibaryons and other strangelet candi- 
dates on the composition and structure of neutron stars including the hyperon degree 
of freedom. We are particularly interested in the question of whether neutron stars 
and their masses can be used to say anything about the existence and properties of 
the H dibaryon. In section |2|, we discuss the condition for the occurrence of dibaryons 
and strangelets in neutron star matter. The relativistic mean field model with hy- 
perons and the H dibaryon is presented in section [3[ Implications for a H dibaryon 
condensate are discussed in section HI and summarized in the last section. 



2 Composite Objects in Neutron Star Matter 

Here we discuss the general features of the appearance of composite quark objects in 
neutron star matter. Nuclei will dissolve in dense matter due to a Mott transition at 
quite low density. Hence, hypernuclei with similar binding energies will also dissolve. 
The situation is different for strangelets, if they are energetically favored compared to 
hadrons. Then strangelets will appear at a certain critical density which will depend 
on the chemical potentials and the mass of the strangelet. The most stable strangelet 
candidates will have a closed shell, i.e. they have zero total spin and are bosons. Also 
the H dibaryon (consisting of two u, d, and s quarks) has zero spin and will form a 
Bose condensate if it appears in dense matter. 

The general condition for a Bose condensation of strangelets is that the effective 
energy must be equal to the chemical potential: 

E* s {k = 0) = m s + U{p) = B ■ fi n - Z ■ fi e (1) 

where B stands for the baryon number and Z for the charge of the strangelet with mass 
m s . The corresponding chemical potentials are \i n = hb and \i e (unit baryon and 
unit negative electric charge, respectively). U(p) is the potential felt by a strangelet 
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in a dense environment. Neglecting interaction and modification of the mass in the 
medium, the threshold condition for the appearance of a strangelet is 



-g- = fl n - Vejj ■ (2) 

Hence, the baryochemical potential determines the onset of condensation as the charge 
to baryon ratio is between +2 and —1 and the electrochemical potential is much 
weaker. 
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Figure 1: The chemical potentials per baryon number for a 
strangelet with a charge of Z/A = —1 (hb + /-O, (/ib), and +1 
(fiB — He)- The horizontal line is the H dibaryon mass of m# = 
2200 MeV assumed to be density independent. The crossing of 
this line with the line for /ib at p c ~ 0.36 fm^ 3 marks the onset of 
H dibaryon condensation. The H mass with a medium dependent 
repulsive potential of Uh = 30 MeV at p is also shown by the 
curve {m H + ?7^)/2. 



We illustrate threshold conditions for the H dibaryon (or any other strangelet with 
assumed free mass per baryon of 1.1 GeV) in Fig. [I] under two circumstances: (1) free 
dibaryon, and (2) dibaryon whose mass is modified by the medium. The two solid 
curves represent the dibaryon mass per baryon under these two circumstances. More 
generally, the chemical potentials for a strangelet with Z/B = —1, 0, and +1 are also 
plotted. The point where these lines intersect and rise above the mass per nucleon 
for a particular strangelet marks the density threshold above which the strangelet in 
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question would comprise one of the constituents of matter. The threshold for free and 
medium modified dibaryons (or strangelets) of free mass per baryon of 1.1 GeV can 
be read from the intersections. Nucleons, leptons as well as hyperons are included 
in the composition of matter. The effect of hyperons is clearly evident in Fig. [I] 
through their saturating effect on p e at p 0.4 fm -3 . The other chemical potentials 
accordingly increase less rapidly with density above the hyperon thresholds. 

The lowest mass strangelet is the H dibaryon which is bound due to color magnetic 
forces |IJ. Assuming the H has a mass of ran — 2.2 GeV and does not change in 
the medium, gives a critical density of about 2p . Note that a negatively charged 
candidate appears at a much lower density and for the assumed mass, would be 
a constituent of matter at densities above p » 0.16 fm~ 3 . If the H dibaryon feels a 
repulsive potential at saturation density of Un{po) = 30 MeV then the critical density 
is shifted beyond the maximum density reached in the interior of a neutron star for 
that equation of state. Note that for a slightly smaller repulsive potential or for 
matter without hyperons the H dibaryon will appear (because its chemical potential 
p,B will rise above its mass as modified by the medium). 

A similar analysis can be done for other strangelet candidates. However, strangelets 
with mass numbers of B < 6 are not stable due to a repulsive color magnetic inter- 
action except for the H dibaryon ||25|| . Nevertheless, strange dibaryon states have 



been predicted to be bound in a relativistic quark potential model fl29| . Negatively 
charged candidates are e.g. and which are heavier than the H dibaryon 

but might appear at a similar density as the H due to their negative charge. It has 
been also proposed in the MIT bag model that negatively charged strangelets with 
closed shells are likely to be most stable against strong and weak emission of hadrons 
P7| . The candidates with a closed shell found are e.g. for B = 10, Z = —4, for 
B — 12, Z — —6, and B — 16, Z — —10. The masses are not precisely known as 
they depend crucially on the value of the bag parameter. To be metastable, their 
masses per baryon should lie between m= and m n . (Absolute stability is unlikely 
because of the finite size shell effect on the quark wave functions.) Mass estimates 
ranges between ms/B = 1.04 — 1.24 GeV for the above candidates. ^From Fig. [I], 
the negatively charged strangelets would appear in neutron stars at densities above 
» 0.4 fm -3 if their masses are lower than 

^ < ^T x + \^T X ~ 1-23 GeV (3) 

(assuming no interactions). Therefore, if strangelets do not feel too high a repulsive 
potential in the medium they can appear as a Bose condensate in neutron stars. In 
the following we will discuss the modification of the properties of neutron stars due 
to the appearance of a strangelet condensate. We choose to study the case of the H 
dibaryon as it is the lightest candidate and might appear first in dense matter. 



3 Mean— Field Model with Dibaryons 

First note that the H dibaryon is a boson with zero spin and isospin. We use the 
standard extended o — uj — p model to describe the baryon sector interacting through 
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the mesons 



B 



2 p 2 



-tf(<r) - \v, u V^ + \mlV»V» + U(V) - + \m%fr (4) 

where B is summed over all states of the baryon octet, the scalar meson is denoted 
by cr, the vector mesons are denoted by and for the iso-scalar and iso-vector 
meson and = dpV v — d u V^. We have taken into account possible self-interaction 
terms for the scalar field U(a) 0] and for the vector field |3l] 

U(a) = l -bm{g a af + \c{g a af , U(V) = -^V^f . (5) 
The H dibaryon is coupled to the mean fields by a minimal coupling scheme following 



20 



C D = V^H*V^H - m* H H*H (6) 
where the vector fields are coupled via the standard replacement 

Vp = d li + ig uH V li . (7) 

This ensures consistency with Ward identities, i.e. the vector fields are coupled to a 
conserved current. The effective mass of the H is defined as in the baryon case 

m* H = m H - g aH a . (8) 

This gives, as for the vector field, a quadratic coupling term of the H to the scalar 
field in the Lagrangian (||). It turns out, that with this choice of coupling, the scalar 
and vector density for the H are the same in the mean field approximation. It was 
also shown that this model is thermodynamically consistent [20]. The equation of 
motion is simply 

D*2r + m* H 2 ] H(oj, k)=0 . (9) 
For s-wave condensation (k = 0) one gets the dispersion relation 

w// = m* H + g uH V = Hh = 2 Hb (10) 

in the mean field approximation, where \xb is the baryochemical potential. This 
relation fixes the amplitude of the H dibaryon condensate. The density of the H 
dibaryon is increased until the effective energy of the H dibaryon is equal to its 
chemical potential. Note that this implies that there must be a repulsive potential 
between the H dibaryons at a certain H dibaryon density. Otherwise, the effective 
energy is decreasing with increasing H dibaryon density and it will never attain its 
chemical potential. 

We do not repeat the full set of equations for the baryons as they can be found 
in e.g. [11J in detail. We note the additional terms due to the H dibaryon condensate 
in the equations of motion 

mla+^-U(a) = Y,9*bP* + 2g aH m* H H*H 
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m 2 . 



V + dV 3 = Y,9»bPv + 19u>h{»h-9»hV q )H*H . (11) 



Here one needs to define only one density for the H dibaryon due to the dispersion 
relation fllPf) 

p H = 2m* H H*H = 2{ t x H -g u}H V )H*H . (12) 
The H dibaryon contributes to the energy density in the form 

e H = 2m* H 2 H*H = m* H p H (13) 

but contributes to the pressure only indirectly through the modification of the meson 
fields via the additional terms in the equations of motion. 

3.1 Baryon-Meson interactions 

There exists various parameterizations in the literature for the nucleon-nucleon inter- 
actions in the mean field model. The parameters are either fixed by nuclear matter 
properties or fitted to properties of spherical nuclei. For example, the parameter 
set used in [|32| with scalar self-interactions U(a) corresponds to the nuclear matter 
properties: B/A = 16.3 MeV, p = 1-53 fm" 3 , a sym = 32.5 MeV, K = 240 MeV, 



m*/m = 0.78 (which for brevity we refer to as GM91). The parameter set TM1 |33 
has been fitted to the binding energy, radii, and surface thickness of heavy nuclei. 
The latter model has a self- interaction term for the vector field U(V). This set has 



been used in p4j| . We adopt these two models as guidelines in the following. 

The hyperon coupling constants have also been chosen differently. We will consider 
two cases: (1) universal coupling of the hyperons and (2) coupling constants using 



SU(6) relations. In the former case, all hyperons are coupled equally |32 



— = ^ = ^ = 0.6 (14) 

9a 9p 9p 

where Y stands for the hyperons A, S, and S. Note that the A has isospin zero, the 
E has isospin 1, while the S and nucleon have isospin 1/2. This gives an additional 
factor of 2 for the p — £ term in the Lagrangian and a vanishing iso-vector coupling 
constant for the A. 



In the other case, SU(6) relations |35| are used for the vector coupling constants 
of the hyperons 

9u ■ 9uA ■ 9u>?, ■ 9cj~ = 3 : 2 : 2 : 1 (15) 

which scale according to the number of light quarks of the baryon. The iso-vector 
coupling constants scale with the isospin like in the universal case but are fixed 
differently to the nucleon iso-vector coupling constant. In our notation this means 

9 P = 9 p t. = 9p~ , 9 p a = . (16) 
The SU(6) symmetry takes already care of the isospin so that the notation as used 

9 P ■ 9 p a ■ 5V£ : 9p~ = 1 : : 2 : 1 (17) 
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means the same. It reflects the strength of the iso-vector potentials of the baryons 
but starts then from a different, SU(3)-symmetric Lagrangian. 

Both cases, universal and SU(6), are consistent with A hypernuclear data insofar as 
the potential depth of the A in normal nuclear matter is fixed to its phenomenological 
value of Ua(po) = —30 MeV. For the universal case the vector coupling constants are 
adjusted to this potential depth for all hyperons. For the SU(6) case the scalar 
coupling constants are adjusted. Note that in addition, the SU(6) coupling scheme 
(quark model) is successful in describing the small A-hypernuclear spin-orbit splitting 

EM- 

3.2 Dibaryon interactions 

The value of the coupling constants of the H to the scalar and the vector field, g a u and 
QuH are unknown. They must satisfy two constraints: (I) the H should not appear in 
normal nuclear matter, and (II) the interaction should allow for neutron star masses 
at least as large as the well established mass of the Hulse- Taylor pulsar. 

The simple quark counting rule suggests that g^n/g^N = 4/3 as the H has four 
light quarks. This choice is motivated from the success of using SU(6) (quark model) 
relations for hyperons in describing hypernuclear properties ||3~5| , [36|. We use it to fix 
g uH . Having fixed the vector coupling without regard to the above constraints (I) 
and (II), the burden of satisfying the constraints falls on the scalar coupling (against 
the background of the other couplings defining the Lagrangian of the theory). 

To fix the scalar coupling g a n we consider the following. A range for the scalar 
coupling constant g a n can be determined from values of the H potential in the medium 
at p 

U H = -gaHcr + guHVo (18) 
because the meson fields are known. The H potential at po should not be deeper than 

Uh(po) > 2E F -m H = 2(m N - 16 MeV) - m H « -350 MeV (19) 

else the H dibaryon would condense at normal nuclear matter density. We choose 
specific discrete values in the range of Uh(po) = +30, 0, -30 MeV. To each, a specific 
value of g a H is implied through Eq. (|TB]). We will find that potentials deeper than 
—30 MeV would decrease the limiting neutron star mass below observed masses (at 
least in the parameterizations of the Lagrangian that we have considered). 

The energy of the H dibaryon is plotted in Fig. [| for neutron star matter including 
hyperons for the various H potentials. The H dibaryon feels a repulsive potential 
above normal nuclear density irrespective of the chosen potential at po- The repulsive 
high density behavior arises from the interaction of the H and the vector meson. 
This repulsion generally dominates at high density unless the scalar interaction is 
too strong. We discuss this issue below. The slope at high density is quite similar 
for the set GM91 and the curves for the different H potentials are just shifted. The 
vector potential dominates at high density and is chosen to be the same (4/3 that 
of the nucleon) in all three cases giving the same slope at high density. The shift 
comes from the differences in the scalar potential which saturates at high density. 
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Figure 2: The energy of the H dibaryon in neutron star matter 
as a function of baryon density with Uh(po) = —30, 0, +30 MeV 
from bottom to top. Solid lines stand for the parameter set GM91 
using universal couplings, dotted for SU(6) couplings, dashed lines 
for the set TM1. 



The behavior at high density is quite different for the parameter set TM1. Here the 
vector potential has a nonlinear dependence on the baryon density due to the vector 
field selfinteraction terms. Note that the H energy and thus the baryon chemical 
potential stays rather constant with density for the lowest curve (Uh(pq) = —30 
MeV). It is especially clear from Fig. |2| that the vacuum mass of the H is not as 
crucial to the appearance of the H in neutron stars as its interactions with vector and 
scalar mesons. 

In place of the above considerations for fixing a range for the scalar coupling 



constant, Faessler et al. 20 invoked the condition 



&h < (2Q) 



ml m? 



on the grounds that otherwise the Yukawa potential between H dibaryons would yield 
a negative compressibility. This is true at low density. For example, for the parameter 
set GM91 we find that H matter is unstable against compression at low density when 
Uh(Po) < +2 MeV. For the parameter set TM1 the low density instability arises 
when Uh{po) < —10 MeV. However, the situation is more complicated. Even if the 
equation of state has a negative slope at low density, it can become positive at higher 
density. Whether or not, cannot be stated in terms of the inequality of Eq. (|20|) but 
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involves all of the other interactions and particle types in matter. Test of the low 
density behavior is therefore insufficient. 

Because gravity compresses a star, the question of stability arises for condensed 
matter, not for low density matter. For the above quoted models the scalar and vector 
potentials have a different (nonlinear) behavior at high density which can alter the 
low-density conclusion. This is demonstrated in Figure [| The binding energy of pure 
H dibaryon matter is shown for the two parameter sets and the three different choices 
of the H dibaryon potential in nuclear matter. One sees that the compressibility, which 
is proportional to the slope of the curves, is negative at low density for Uh(po) — — 30 
MeV for both parameter sets. Nevertheless, the scalar self-interactions provide a 
nonlinear dependence of the scalar potential on the density (here H dibaryon number 
density). This results in an overall repulsive potential for pure H matter at higher 
density and a minimum around normal nuclear density. Hence H matter can be stable 
at high density even if the low density limit seems to indicate an instability. 




Density [fm~ 3 ] 

Figure 3: The equation of state of pure H matter for various 
potential depths in nuclear matter (Uh(po) = —30,0, +30 MeV 
from bottom to top). Solid lines stand for the parameter set 
GM91, dashed lines for TM1. 

With reference to Fig. |3], we see that even more complicated situations can arise. 
For the parameter set TM1, the vector self-interactions cause the energy density due 
to the vector potential to rise like p 1//3 instead of like p, as for the linear behavior 
in the standard Walecka model. It can then happen that the scalar potential wins 
ultimately over the vector potential at very high density. Such a case is seen in 
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Figure |3] for the lowest dashed curve: the equation of state has a local minimum 
around normal nuclear density but the compressibility becomes negative again at 
higher density. 

The instability for pure H matter does not mean that a neutron star with a H 
condensate is unstable against compressional modes, as the overall compressibility can 
still be positive. This will depend on the intrinsic stiffness of the equation of state 
used and the hadron population inside the neutron star. A recent analysis considered 
neutron star matter consisting of nucleons and leptons only [[H]] and found rather 
stringent conditions for the coupling constant of the H dibaryon. The appearance of 
hyperons already at (2 — 3)p HU will certainly alter their conclusions and will be 
discussed in the following. 



4 H Dibaryons in Neutron Stars 

4.1 Populations 

To give an early impression of the possible presence of H dibaryons in neutron stars 
under acceptable conditions as to (1) its absence in normal matter, and (2) an ac- 
ceptable value for the limiting mass neutron star, we compare the populations of the 
limiting mass star in the absence of the dibaryons and in their presence, in Figs. ^ and 
H Hyperons appear abundantly in the interior of the neutron star. The hyperons A 
and S~ reach values close to the proton density in the stellar core. Protons are more 
abundant once the negatively charged hyperons S~ and S~ are present to compensate 
the charge of the proton. As can be seen in Fig. ||, particle populations are strongly 
modified in the core of the star where the dibaryon appears. The proton population 
is suppressed since baryon number is carried more in the H bosons. Likewise the 
hyperon populations are strongly suppressed. This is not to say that these baryons 
have little influence on the H presence. It is in matter containing significant S~ and 
A populations that the H threshold is attained. Outside the core, beyond r w 4 km, 
populations and the stellar radius are hardly changed by the presence of the H in the 
core and therefore cooling of the star would be little effected by the H. H dibaryons 
could lie within the star while providing no direct sign of their presence. 

4.2 Limits on the H dibaryon mass 

Generally the maximum mass of a neutron star is lowered due to H dibaryon con- 
densation because the condensate does not contribute directly to the pressure and it 
removes the pressure due to two baryons per dibaryon in the condensate. The equa- 
tion of state is thereby softened. As the maximum mass should be not lower than 
1.44M (the mass of the Hulse- Taylor pulsar), one can try to impose constraints on 
the medium mass of the H dibaryon. For the nuclear parameterization TM1, we 
show the equation of state corresponding to several values of Uh(po) in Fig. Q The 
equation of state can be considerably softened by the condensation of H dibaryons. 
Especially for the case Uh(po) = —30 MeV the pressure stays nearly constant once 
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Figure 4: Populations of octet 
baryons and leptons in a limiting mass 
star with nuclear properties as de- 
scribed by the case GM91 with hy- 
peron couplings chosen as in the SU(6) 
scheme. 



Figure 5: Similar to Fig. ^ but with 
the H dibaryon experiencing a poten- 
tial of Uh(po) = -30 MeV. 



condensation starts. The H boson constitutes a large fraction of the matter but does 
not contribute to the pressure directly. The plateau seen in the equation of state in 
Fig. |6] can be traced back to the equation of state for pure H dibaryon matter in Fig. |3| 
(see the bottom dashed curve). 

Figure [7] summarizes the neutron star masses for the two hyperon coupling schemes 
SU(6) and 'universal' (see section |3~T|) and for various values of the interaction Uh(po)- 
The more attractive the interaction, the more populous the dibaryon, the softer the 
equation of state and the smaller the limiting mass. The potential Un{po) = —30 
MeV is about as attractive as is compatible with the Hulse- Taylor pulsar. For the 
SU(6) coupling, H dibaryons feeling Uh(po) = MeV would not be present in the 
stable members of the sequence, while for universal coupling of hyperons, the dibaryon 
is present in small number and reduces the limiting mass marginally by ~ 0.03M Q . 
A repulsive dibaryon interaction in the medium would therefore assure its minimal 
presence if not its total absence (as is the case for the SU(6) coupling). 

One might infer from the figure, that very attractive potentials for the H of 
Uh(po) < —30 MeV can be ruled out by neutron star data. Nevertheless, such a 
conclusion must be moderated by our ignorance of the H dibaryon mass. The ap- 
pearance of the H condensation depends on the in-medium potential at high density 
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Energy Density [MeV fm~ ] 

Figure 6: Equation of state for the parameter set TM1, SU(6) 
coupling for hyperons and several values of the H dibaryon po- 
tential Uh(po)- 



and the mass of the H dibaryon which are unknown. 

In addition to uncertainties in the vacuum and medium mass of the H, there are 
uncertainties in the underlying nuclear equation of state especially at densities above 
normal nuclear density. This can be seen by contrasting Fig. |7| with Fig. [8[ In the 
second case the nuclear parameterization is TM1 of section [3.11 Here we see that the H 
is present provided Un(po) < +30 MeV. For an attractive potential of Uh(po) = —30 
MeV, one finds that the mass of the neutron star reaches a plateau, i.e. the mass of 
the neutron star is independent of the central energy density. The maximum mass is 
then 1.44M . This is the lowest value allowed by present observation. The plateau 
comes from the fact, that the equation of state has a nearly constant pressure for 
the particular coupling due to the appearance of the H dibaryon. A more attractive 
potential than Uh(pq) = —30 MeV would result in a negative curvature of the pressure 
and therefore in an unstable equation of state. Note, that this behavior is mainly 
related to the potential between the H dibaryons as discussed in connection with pure 
H dibaryon matter (see Figure ^). Neutron star matter is stabilized against collapse 
to H dibaryon condensate because of the repulsive interaction with the vector meson 
provided the H dibaryon density is not too large. The presence of the other baryons 
in the matter tends to stabilize the system compared to pure H matter. 

In the present approach one can now exclude certain regions of the assumed mass 
of the H Dibaryon mj? and its potential at saturation density Uh- For a too low H 
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Figure 7: Details of neutron star se- 
quences near the limiting mass for the 
nuclear parameterization GM91 and 
two hyperon coupling schemes labeled 
'universal' and SU(6) and for each of 
these, several values of the dibaryon 
interaction Uh{po)- 
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Figure 8: Similar to Fig. [7| but for the 
nuclear coupling denoted by TM1. 



mass or a too deep potential the presence of the H dibaryon condensate will lower the 
maximum mass of a neutron star below the observed limit of 1.44M . This excluded 
region is shown as a shaded area in Fig. || in an (£/#, run) diagram for the parameter 
set GM91 using universal coupling. The dashed line denotes the case for the parameter 
set TM1. The thresholds for AA and EN decay which are relevant for the lifetime of 
the H dibaryon |37J are also indicated. The results for the two parameterizations are 
quite close to each other despite their different high density behavior. For an H mass 
of lower than 2.13 GeV (the EN threshold) an attractive potential in nuclear matter 
gives too low maximum neutron star mass. Hence, rather long-lived H dibaryons 
(r > 10~ 7 s according to |37|]) are unlikely to form bound H hypernuclear states. On 
the other hand, attractive potentials lower than Uh{po) < — 50 MeV seem to be ruled 
out by neutron star mass constraints for the H mass range considered. Otherwise, 
the H dibaryon has to be heavier than the 2A threshold and will be a resonance state. 
These limits will depend also on the chosen vector coupling constant. We are using an 
effective model and the extrapolation to high density might be completely different 
in reality. 
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Figure 9: Diagram of the H vacuum mass and the H binding 
energy in normal nuclear matter. The shaded region gives a max- 
imum neutron star mass lower than 1.44M Q for the parameter set 
GM91 using universal coupling, the dotted line denotes the case 
for TM1. 



4.3 Radius and the Mass-Radius Relation 

The mass-radius relation for the GM91 model with universal hyperon couplings is 
shown in Fig. |l(] for several values of the dibaryon potential Uh(po)- The more at- 
tractive the potential, the softer the equation of state and the lower the limiting mass, 
as remarked earlier. The radius of the limiting star decreases the more attractive the 
potential Uh(po) because the star, having less mass, is less gravitationally compacted. 
Similar results are shown in Fig. [11] for the nuclear parameterization TM1. In this 
case the limiting mass stars have substantially larger radii. 

It is interesting to note, that the presence of the H dibaryon in neutron stars seems 
to lower the maximum mass but increases the minimum radius. The mass-radius 
relation just stops at the point where the H dibaryon condensation sets in. This is 
contrary to kaon condensation |38|, |39| where the radius decreases for a kaon condensed 
star. Note that kaon condensation produces a strong phase transition of first order but 
the equation of state has no plateau if treated in a thermodynamically consistent way 
40| . In addition, baryons are not replaced by a baryon number carrying condensate 



but neutrons are replaced by protons and K 's resulting in a different mass-radius 
relation. 
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Figure 10: Mass-radius relation for 
the nuclear model GM91, universal hy- 
peron coupling and several values of 
the dibaryon potential Uh(po)- 
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Figure 11: Similar to Fig. [L(J but for 
the nuclear model TM1. 



5 Summary 

We are particularly interested in the question of whether neutron stars and their 
masses can be used to say anything about the existence and properties of the H 
dibaryon. We have studied the influence of the possible occurrence of an H dibaryon 
condensate and strangelets in neutron stars including hyperons. Without in-medium 
modifications, it is quite likely that especially negatively charged strangelets, if they 
exists, will be present in the dense interior of neutron stars. 

The appearance of H dibaryons in the stellar core depends crucially on their mass 
and on the chosen potential of the H in nuclear matter. Hyperons tend to shift 
the onset of the H to higher density or to prevent H dibaryon condensation. If the 
condensation happens and if the potential of the H is attractive enough to provide a 
substantial number density in the neutron star, the maximum mass of the neutron 
star is reduced compared to the case without the H dibaryon. The decrease of the 
maximum mass is moderate and allows for the presence of H dibaryons in the interior 
of neutron stars in accord with present neutron star mass data. If the H dibaryon 
feels an attractive potential in matter, it can lead to a plateau in the mass of neutron 
stars, as there exist a region of very slowly rising pressure with energy density. 

If the limiting neutron star mass is about that of the Hulse- Taylor pulsar a conden- 
sate of H dibaryons of vacuum mass ~ 2.2 GeV and a moderately attractive potential 
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in the medium could not be ruled out. On the other hand, if the medium potential 
were even moderately repulsive, the H, would not likely exist in neutron stars. If 
neutron stars of mass ~ 1.6M Q were known to exist, attractive medium effects could 
be ruled out. For a mass limit of 1.44M Q , attractive potentials for an H mass below 
the EN threshold (1.3 GeV) are ruled out. 

H dibaryon or strangelet condensation might happen as a precursor to the phase 
transition to a quark plasma. In this respect, we note that this phase transition is 
of first order [[|lj. Hence, small bubbles of strange quark matter will appear in the 



mixed phase which are most likely negatively charged due to the isospin potential 
of the nuclear matter. This is in line with the results presented here. As the most 
stable strangelets have spin zero, the onset to a quark plasma will be initiated by a 
Bose condensation of strangelets (possibly including the H dibaryon). As the phase 
transition proceeds, the bubbles will overlap and will finally replace nuclear matter 
by essentially filling up the whole volume. 
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